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Abstract. A new pairing theory for many-fermion systems is obtained via the Dirac 
super symmetry framework recently introduced to describe Dirac particles in external 
potentials. It is shown that the standard Bogoliubov-Valatin canonical transformation 
treatment of the quasi-particle BCS singlet pairing mechanism naturally falls within 
this framework. Straightforward generalizations in which the fermions can be ascribed 
v components are shown to lead to enhanced gap energies and critical temperatures 
as in the case of cuprate superconductors without invoking a stronger electron-boson 
coupling. The new T™ ax limit is T™ ax = vTf cs ', with T? cs » 40° K. 
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One of the most striking approaches to the problem of pairing was pioneered by 
Nambu |l| who showed that a remarkable connection exists between the pairing 
mechanism and supersymmetry ||. Within this framework, supersymmetry must be 
severely broken in order to make contact with phenomenology. It has recently been 
shown, see Refs. Jl, £|, [|, that a generalization of quantum supersymmetry to massive 
Dirac fermions || can be achieved if the Dirac Hamiltonian is of the form 

H = Q + Q^+A, (1) 

with A a Hermitean operator, provided that these operators satisfy the following 
anticommutation relations 

{g,A} = {g t ,A} = o = g 2 . (2) 

This new structure has been named Dirac supersymmetry (Dirac Susy). 

In this letter we apply Dirac Susy concepts to superconductivity phenomena 
and demonstrate that one may accommodate within this formalism both the usual 
BCS mechanism as well as other kinds of fermion pairing. In view of the generality 
of the approach, it is likely that the underlying mechanisms of the so-called exotic 
superconductors can be treated within the scope of Dirac Susy. One indication of 
this is seen in the example of a two-layered fermion gas in which a substantial increase 
in the energy gap can be achieved even for weak interaction. Our formalism lends 
itself to study the consequences of a BCS mechanism for a general fermion system 
with several components within the canonical transformation scheme proposed by 
Bogoliubov and Valatin (BV) [§]. For a recent survey of the relevance of BCS-based 
mechanisms to high-temperature superconductivity (HTSC) see Ref.f9j. The notion of 
component that we use should be understood in a broad sense: it can arise physically 
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from a spin degree of freedom, or from a pseudo-spin generated by layers in which an 
electron or hole moves; it can also stem from a finite number of transversal fermion 
states in a slab of finite thickness; or in general from any condition leading to a 
discrete representation of an degree of freedom. The latter is the case of the electron 
gas [0 in a doped multi-valley semiconductor or multi-valley semimetal. 

A Dirac Susy structure such as ([1]) admits an exact unitary Foldy-Wouthuysen 
(FW) transformation of the original Hamiltonian || ^, |5| . The squared Hamiltonian 
becomes if 2 = {Q, Q^} + A 2 = h 2 + A 2 , where h is required to be an even root of the 
supersymmetric operator, h 2 = {QjQ^} in the FW sense. Defining A = A/v^A 2 , the 
FW transformed Hamiltonian is then 
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H FW = A^{Q,Q^} + A 2 , (3) 
provided A has non-null eigenvalues. Further, Q and in matrix form are 

(4) 

where the q are in general submatrices. Because A commutes with Hfw an d its eigen- 
values are ±1, Eq. @ implies generally that, the energy spectra has two branches; 
between them there is a gap determined by the lowest eigenvalues of A 2 and h 2 . If 
A has negative eigenvalues and H 2 is not bounded from above one gets an unstable 
ground state for the Dirac Susy Hamiltonian (]I|). The usual remedy is to introduce a 
Dirac, or Fermi, sea. The important feature of Dirac Susy interactions is that the def- 
inition of this sea is coupling strength independent, this property being the stability 
of the Dirac sea || . The introduction of the sea necessarily implies a field theory. 

Note that a standard quantum mechanical Susy is approximately obtained if A 2 
eigenvalues are either very large or very small as compared to those of h 2 . Therefore, 
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Dirac Susy is a specific type of supersymmetry-breaking. It is more appropriate 
to consider it as a generalization of the usual supersymmetry to cases when the 
Hamiltonian is fermion-like. A relevant consequence to our purpose here is that 
Dirac Susy restores and guarantees a gap in the energy spectrum of the system. 

Let us first consider how this formalism is connected to the BV transformation 
in BCS pairing phenomena. The essential characteristic of BV theory is to exhibit the 
well-known Fermi sea instability \\ll\\ through a canonical transformation which mixes 
fermion (electron, for definiteness) and hole states of opposite momentum and spin 
quantum numbers ( viz., k j and — k !). The physical basis for this transformation 
relies on the fact that individual electrons in the Fermi sea near the Fermi surface 
escape the sea and become bound into Cooper pairs. The new states, so-called bo- 
golons, are related to the electron states, and a_kj,, by a transformation matrix 
which must be orthogonal in order to preserve the anticommutation relations for the 
bogolon operators, «k and /3_k- We follow the notation of Fetter and Walecka ||12|| . 
The new vacuum state |0> satisfies «k|0> = /5k|0> = 0. A realization of this new 
vacuum could be, e.g., the BCS state |0> = M fli (1 + 9i a i a -i)\ ( t ) >> where \<p > is 
the original (quasifermion) vacuum, is a normalization factor and is related to 
the amplitude of pairing in the % and — % states. Here, i stands for the relevant single- 
particle quantum numbers. Clearly, the new vacuum, |0>, is not annihilated by the 
quasifermion operators a,. At zero temperature the thermodynamic potential ]I2| 
Q(T = 0, V, n) is given by the expectation value of the operator K = H — fiN 

K = H( e ° - ^) a kA a kA (5) 

kA 
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-- kiAik 2 A 2 |l / |k3A3k4A4 > aj CiAi aj C2A2 a k4 A 4 ak 3 A 3 - 
ki,k 2 ,k 3 ,k 4 

Ai, A2, A3, A4 

Applying the BV transformation to this leads to separation of the operator K into 
a zero-body, H , a one-body, Hj, and a two-body, Hu, operators; the weak-coupling 
assumption allows treating Hu as a small perturbation to be neglected. Hi written 
in a Nambujl]] matrix form is 
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(6) 
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a^ — ik 

where measures the Hartree-Fock quasi-particle energy e° with respect to the chem- 
ical potential /1 and A& is the gap function A^ = J2k' < k — k\V\k' — k' > UkVk- As 
will be shown shortly, @ already implies a Dirac Susy structure. 

The BV transformation is fixed by the requirement that this one-particle Hamil- 
tonian Hi should be diagonal for a nonzero value of Vk with respect to the a, (3 basis. 
The bogolon excitation energy is then given by Ek 



£| + A|, and the gap function 
Ak is determined from the self-consistent equation A& = | J2k' < k — k\ V\ k' — k! > A&/ 
The superconducting solution being characterized by A& 7^ 0. The familiar BCS 
interaction model giving a simple yet realistic nontrivial solution is of the form 
< k - k\V\l - I > = (V /L d )9(hu D - \ik\)6(hu D - |^|), where V is the (square 
of the) electron phonon coupling, L d in d dimensions, is a normalization volume and 
ujd is the cutoff (Debye) frequency . Introducing iV(0), the density of states for one 
spin projection at the Fermi surface, one obtains the celebrated gap equation, 
V o N(0) [ h "D d£ 



V N(0) sinh-^^/A), 



(7) 
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Letting VoiV(O) = g and solving for A one gets A = hcu D / sinh(l/p) ~ 2huc,e^ 1 ^ 9 , 
for weak coupling g <C 1. For pure elements g is in the interval 0.15 < g < 0.6. In a 
finite-temperature formalism this implies T c w 1.13fuvD^~^ • 

We now generalize the BV transformation for a system in which each fermion is 
characterized by a discrete parameter A that can take on v values. For v = 2, this 
could correspond to the two spin-| degrees of freedom. Alternatively, if the electrons 
are confined to a slab of thickness t, the discreteness in A corresponds to the different 
excitations of the transverse degree of freedom and in addition to the spin. Other 
concrete examples are a many layered electron gas, for which A refers to the specific 
layer, and the multi-valley structures referred to before. In close analogy with 
the one-component case let us now define 



{ a kl, 



O'kui a -k\i ' ' ' O'-k 



(8) 



where T stands for transpose matrix, and 



B, 



\ 



A k = UA k 



(9) 



where vxv matrices and //-vectors are denoted with a check above. The commutation 
relations for the operators ak\ can be reexpressed in matrix form as 



A k Al,±A* k ,Al = 16 kkl , 



(10) 



while the canonicity condition of the BV transformation requires that BkB\,±Bl,Bl = 
l^fcfc'; which is satisfied if U in (||) is an orthogonal matrix, as can easily be verified 
by sandwiching (10) between U and W. When this transformation is applied to the 
Hamiltonian K, one obtains as before an operator of the form K = Ho + Hj + Hu. 
Again the term Hjj is neglected. One has for Hj the form (0) but with all quantities 
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bearing a check, i.e., they are v X v matrices. For our purposes the relevant structure 
is just 



£fc A fe 



If in ([!]) we identify 

Q 



q 



/ 0^ 



-6 



A fc 



and A 







(12) 



/ 



one can easily confirm our central result that the standard BV operator and its extrap- 
olations to multicomponent systems have in fact a Dirac Susy structure, and therefore 
it guarantees an energy gap. In this case, the BV transformation is completely equiv- 
alent to the FW transformation. Quite generally the condition (Q) implies 



[A,|] = o, 



(13) 



which is satisfied if £ = £q1 + £i A , where £o and £i are functions of k. 

Let us now introduce a finite-temperature formalism, it is useful to consider 
a mean-field approximation |T3| , [14|]. The thermal expectation values are < A >= 
Tr(e~ f3K A)/Z, with Z = Tr[e~^ K \ one has in this context that < ajaLj > ^ ^ 
< did-i > where we have merged the indices k and A into the single i. With this 
notation, the matrix elements of A in the finite temperature formalism are 



AaiA 2 (&) = < ^ ^2 - k\V\X 3 l; A 4 - / >< a_ lXi ai Xs > 
An effective one-body Hamilton operator can be written as 



(14) 



K = [e(/c) — fi]a\ai 



(15) 
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k 



A k , (16) 



which leads to the Dirac Susy structure in the finite temperature formalism 

[e(k)-/i]l A(k) ■ 

as in the Hi case. Then we have found Dirac Susy structure in the finite temperature 
formalism in a system with multicomponentes, independent of the structure and the 
intensity of the coupling. 

We now illustrate the general ideas developed above with the following examples 
with v = 2 for which the components are not (physically) equivalent to spin. This in 
practice means that the spin degree of freedom will be assumed to couple to a singlet 
state in the Cooper pair. The physical systems to which these models correspond are 
in fact inspired in the quasi-2D electron gas. a) First, consider an electron gas confined 
to a 2D slab of finite thickness; we assume that only the two lowest transversal energy 
states are important, i.e., the energy separation between the successive transversal 
states is big enough so that higher states can be neglected, b) Secondly, consider a 
two-layer 2D electron gas for which the spin degree of freedom could be neglected, and 
take the case for which the fermion degrees of freedom have been diagonalized before 
the BV transformation. These therefore describe quasiparticles a la Landau, and as 
emphasized by Schrieffer [[0]] the pairing mechanism is more adequately applied to 



these states. Formally these states are described by diagonal matrices In both 
these models we take = l£fc! physically this means that the two components of 
(quasi-) electrons are equivalent before the BV transformation. Consequently ( pT3|) is 
satisfied for an arbitrary matrix A*.. Performing the BV transformation on ( |T5D we 
derive the bogolon Hamiltonian 
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This form is further diagonalized by any transformation diagonalizing A^A^ and 
Aj,Afc, it being sufficient to study either one of these two matrices because they 
appear in diagonal blocks and have the same eigenvalues. It is convenient to expand 
A in terms of Pauli matrices, f , namely: A = A 1 + A • f , with A a 3- vector, where 
the subindex k has been omitted, and Ao and A might be complex. Then, 

AA f = [A A* + A • A*] 1 

+ [(A A* + A*A)+*A x A*]-f, (18) 

where the term proportional to f is also real. Furthermore, the two vectors (A A* + 
AqA) and A x A* are clearly orthogonal to each other. Consequently the diagonal 
form of AA^, (AA')d, is simply 

(AA% = [A A* + A ■ A*]l + ^(A A* + ASA) 2 + |Ax A*| 2 f 3 . (19) 

Now, the lowest (highest) transition temperature T c is essentially determined by the 
smallest (largest) eigenvalue of A*., via A(T C ) = 0. Eq.(|19D is the most general form 
of the gap equation for the effective Hamiltonian of Eq . (|i~6l) , with v = 2. 

Let us now discuss some concrete realizations of this formalism. First, it follows 
that the conventional BCS result can be obtained with A = and A = 5. This gap 
form is consistent with Vx^XsX^ki = 5\ x x 2 8x s x^/u- For this case a finite temperature 
formalism gives the usual phase diagram and the ratio R = 2A(T = 0)/(kbT c ) = 3.1, 

A different Ansatz which is important in connection with the exotic heavy 
fermion phenomena Q] is found assuming Ao = and real A (up to an overall 
phase). Such a situation can be obtained dynamically with an electron-electron inter- 
action of the form V^a^a^z = 1/2(^3^4 ~ $\ 1 \ 4 5 X2 \ 3 )V kl . If all the components 
are of the same size 5 (which physically means that all fermion couplings leading 
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to Cooper pairing are equivalent) we get an enhancement of the gap for this triplet 
interaction, A(0) = \f?>5. In this case one gets a phase space diagram which is of the 
same generic form of the traditional BCS but for which A(0) = \f3A(0)Bcs', that 
implies an increased T c = \/3T c bcs- A new physical situation arises for the case of 
Ao = Aj but for this the enhancement is greater, T c = 2T c bcs- This means that 
triplet and singlet Cooper pairs are combined and that they add constructively to 
T c . With a greater number of components, v, and assuming a equal coupling among 
them, we conjecture that even larger gap enhancements of order v can be obtained if 
a material in which all the components interact with similar couplings. 

A different effect arises if one seeks a big splitting in (|L9|) . This is exemplified 
with Ao = and one cartesian component of A being imaginary relative to the other 
two, for example Ai = A2 = 2A3 = 5, with 5 real. Under such conditions there 
are two well-differentiated values for the gap (3 ± \^8)5 2 ; which in turn imply two 
transition temperatures. In a model with more components one expects a multi-gap 
system with multi-transition temperatures. In this discussion we have mainly been 
concerned with discrete fermion degrees of freedom. We remark, however, that this 
generally has implications for the orbital pair wave functions. For example, in the 
layered realization mentioned above, coupling among fermions in different layers will 
require a spherically asymmetric wave function, which in turn will imply P, D, etc. 
wave components p~5f] . 

In conclusion, the general Dirac Susy is adequate for pairing phenomena. There 
exists a general connection between Dirac supersymmetry and the BCS theory of 
superconductivity. Dirac supersymmetry was sketched, in particular the remarkable 
result that a unitary FW transformation decoupling positive and negative energy 
states can be explicitly constructed. The relation to conventional non-Dirac super- 
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symmetry, with a possible central charge extension, was elucidated. It is implicitly 
shown that a much more complex supersymmetry breaking term is needed to under- 
stand superconductivity than previously thought. However, this is naturally taken 
into account by Dirac supersymmetry. So for the case of a singlet conventional pair- 
ing theory we showed that the effective quasi-electron Hamiltonian is precisely of the 
Dirac supersymmetric form, and that the BV transformation lends itself naturally 
to its derivation. The relation between superconductivity and Dirac supersymmetry 
was then generalized to a multicomponent fermionic system. For a two-component 
system it was shown that a simple Ansatz leads to an increase in T c by a factor of 2, 
and for a v component system the gap increase can conjectured to be of order v. 
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